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Abstract. In this study, we evaluate the ability of the 
latest NASA GISS composition-climate model, GISS-E2- 
PUCCINI, to simulate the spatial distribution of snow BC 
(sBC) in the Arctic relative to present-day observations. Ra- 
diative forcing due to BC deposition onto Arctic snow and 
sea ice is also estimated. Two sets of model simulations are 
analyzed, where meteorology is linearly relaxed towards Na- 
tional Centers for Environmental Prediction (NCEP) and to- 
wards NASA Modern Era Reanalysis for Research and Ap- 
plications (MERRA) reanalyses. Results indicate that the 
modeled concentrations of sBC are comparable with present- 
day observations in and around the Arctic Ocean, except for 
apparent underestimation at a few sites in the Russian Arc- 
tic. That said, the model has some biases in its simulated 
spatial distribution of BC deposition to the Arctic. The sim- 
ulations from the two model runs are roughly equal, indi- 
cating that discrepancies between model and observations 
come from other sources. Underestimation of biomass burn- 
ing emissions in Northern Eurasia may be the main cause of 
the low biases in the Russian Arctic. Comparisons of mod- 
eled aerosol BC (aBC) with long-term surface observations 
at Barrow, Alert, Zeppelin and Nord stations show signifi- 
cant underestimation in winter and spring concentrations in 
the Arctic (most significant in Alaska), although the simu- 


lated seasonality of aBC has been greatly improved relative 
to earlier model versions. This is consistent with simulated 
biases in vertical profiles of aBC, with underestimation in the 
lower and middle troposphere but overestimation in the upper 
troposphere and lower stratosphere, suggesting that the wet 
removal processes in the current model may be too weak or 
that vertical transport is too rapid, although the simulated BC 
lifetime seems reasonable. The combination of observations 
and modeling provides a comprehensive distribution of sBC 
over the Arctic. On the basis of this distribution, we estimate 
the decrease in snow and sea ice albedo and the resulting 
radiative forcing. We suggest that the albedo reduction due 
to BC deposition presents significant space-time variations, 
with highest mean reductions of 1.25 % in the Russian Arc- 
tic, which are much larger than those in other Arctic regions 
(0.39 % to 0.64 %). The averaged value over the Arctic north 
of 66° N is 0.4-0.6% during spring, leading to regional sur- 
face radiative forcings of 0.7, 1.1 and l.OWm" 2 in spring 
2007, 2008 and 2009, respectively. 
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1 Introduction 

Emissions of black carbon (BC) particles result from incom- 
plete combustion during the burning of biomass and fos- 
sil fuels, and are considered a significant climate forcing 
factor (IPCC, 2007; McConnell et al., 2007). In the atmo- 
sphere, the absorption of sunlight by BC contributes to global 
warming and alters cloud-formation processes (Jacobson, 
2001). After deposition onto snow and ice, BC has the po- 
tential to significantly reduce the surface albedo, hence per- 
turbing the radiative balance and possibly leading to earlier 
snowmelt (Hegg et al., 2009). The global average radiative 
forcing from BC by altering surface albedo was estimated 
as +0.1 W m -2 (IPCC, 2007), with estimates varying from 
0.01 to 0. 16 W m -2 (Flanner et al., 2007; Hansen et al., 2004; 
2007; Koch et al., 2009a). 

Arctic climate is especially vulnerable to BC deposition 
because of the abundant and wide distribution of relatively 
clean sea ice, snow and glaciers, which can be impacted ef- 
fectively by accelerating melting and snow/ice albedo fee- 
backs (McConnell et al., 2007; Koch and Hansen, 2005). Re- 
cent research suggests that the seasonally averaged surface 
forcing by BC in the Arctic can be up to +0.5 W m -2 in 
spring (Flanner et al., 2007; Quinn et al., 2008). The com- 
parisons between modeled and observed aerosol BC (aBC) 
at Barrow, Alert and Zeppelin stations have shown that most 
previous models underestimated concentrations of BC in the 
Arctic, especially in winter and spring (Shindell et al., 2008; 
Huang et al., 2010a; Liu et al., 2011). However, the sea- 
sonal cycle and magnitude of other trace species, like CO and 
ozone, are reproduced much better in most models (Shindell 
et al., 2008). 

There have been many studies on the dominant causes of 
the present discrepancy between models and observations. 
Several factors were suggested, such as the uncertainty of the 
emission inventory used in the model (Wang et al., 201 1), the 
modeling of transport processes (Liu et al., 2011), the lack 
of significant number of observation sites available to ade- 
quately validate models, or the existing challenges in quanti- 
fying BC mass concentrations by current measurement meth- 
ods introducing an inherent uncertainty on reported BC mass 
concentrations (Shindell et al., 2008; Doherty et al., 2010; 
Baumgardner et al., 2012). In recent years, model simu- 
lations and measurement methods have been updated and 
greatly improved, so that it is useful to collate comprehen- 
sive pan-Arctic sBC observations from previous campaigns 
to re-evaluate the current models. Koch et al. (2009b) com- 
pared vertical profiles of aBC from several models to aircraft 
observations during the International Polar Year (IPY), and 
suggested that current models underestimated BC concen- 
trations throughout much of the troposphere in different de- 
grees. There still exist large uncertainties in this comparison, 
because the observations of vertical profiles of aBC present 
just a “snapshot” of the BC distribution in limited sites and 
they are not representative over the whole Arctic, especially 


at climate scales. More extensive and long-term observations 
are still needed to give a comprehensive validation to cur- 
rent models. This study focuses on the latest NASA GISS 
composition-climate model (hereafter GISS-E2-PUCCINI), 
aiming to investigate the performance of the simulation of 
spatial distribution of BC deposition on Arctic snow and ice. 
We summarize BC measurements reported to date and add 
the observations obtained in the first Korean Arctic cruise 
(2010) to the analysis. We use these to evaluate the modeled 
distribution of sBC in the Arctic and then give a comprehen- 
sive map of the spring sBC in and around the Arctic Ocean 
with a combination of observations and validated model re- 
sults. Finally, we estimate the decrease in snow and sea ice 
albedo and resulting radiative forcing due to BC deposition 
based on the derived sBC distribution. 


2 Observations 

The earliest observations of sBC started from the mid- 
1980s (Clarke and Noone, 1985), and were carried out 
at only a few sites, such as Camp Century, Greenland 
(77.2° N, 61.1° W), Dye 3, Greenland (65.2° N, 43.8° W), 
Alert, Canada (83.5°N, 62.5° W), Barrow, Alaska (71.3° N, 
156.6° W) and several Arctic Ocean sites. Later, the spa- 
tial distribution of BC in snow and sea ice was investigated 
during the SHEBA experiment (Grenfell et al., 2002). Dur- 
ing recent years, the circumpolar regions were surveyed as 
a whole, expanding and updating the previous observations 
(Doherty et al., 2010; Forsstrom et al., 2009), which includes 
near one hundred sites and two thousand samples in the Arc- 
tic and sub-Arctic regions. 

We summarize the spring measurements of present-day 
(mainly from 2007 to 2009) sBC in the Arctic from all stud- 
ies known to us. Some of these include vertically-integrated 
measurements or measurements at various depths throughout 
the snowpack, and others just include surface measurements. 
Concentrations supplied by Doherty and others have been in- 
creased by 1 1 % relative to their earlier study (Doherty et al., 
2010), in order to correct for a low bias in the concentration- 
dependent parameters ( CJj' ( ax , Cg£ and C^ 1 "'' ) which resulted 
from an error in the earlier data analysis. As reported by Do- 
herty et al. (2010), the vertical distribution of sBC in the Rus- 
sian Arctic is non-uniform, with larger values in the upper 
25 % of the snowpack at most measurement sites (left panel 
in Fig. 1). In contrast, the concentrations observed in the 
Canadian Arctic are smaller and relatively uniform through- 
out the snowpack (right panel in Fig. 1). Both sets of vertical 
profiles involve measurements of sBC and snow depth from 
dozens of snow pits, thus they are considered to be largely 
representative of the typical distribution of spring sBC in 
these regions. The typical values of snowpack depth in Rus- 
sia and Canada were respectively 34 + 7 cm and 30 + 14 cm 
across all measurement sites. We calculate the surface and 
subsurface concentrations of sBC from the observations at 
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Fig. 1 . Vertical profile of sBC concentrations in the Arctic snow from the field measurements in west (2007) and east (2008) Russia (left), 
and profiles from the field measurements in Canada in spring 2009 (right). Cg^, denotes the estimated true mass of BC per mass of snow. 
These two figures are respectively from Fig. 12a and Fig. 9 (left panel) in Doherty et al. (2010), but with some necessary modifications to 
highlight different concentrations of sBC in various depths. The snowpack depths in Russia and Canada were respectively 34 + 7 cm and 
30+ 14 cm across all measurement sites. 


different depths, and then give an equation for estimating the 
integrated-layer concentration (C/,, h — hi + h2) of sBC. In 
the case of surface sBC > subsurface sBC (for most measure- 
ment sites in Russian Arctic): 

_ 25%p/,iC;,i +15%piaCi,2 
h ~ 25% phi +75%p; i2 

where Chi is the concentration of sBC in the surface snow 
(h i —25% -h) at each site, phi is the average density of sur- 
face snow, C /,2 is the concentration of sBC in subsurface 
(/?2 = 15% ■ h ), pi ,2 is the average density of subsurface snow 
and hi is the top 25 % of the snow pack depth and h 2 is the 
bottom 75 %. Note that it is difficult to calculate the precise 
Ch based on present limited observations, because we have 
not obtained the snow density values that correspond to BC 
measurements in various snow depths. In this study, we ap- 
ply the estimated snow density of the surface and subsur- 
face layer: p/,i = 0.256gcm -3 and p /,2 = 0.345 gcm~ 3 that 
is calculated from the mean values of snow density in differ- 
ent types of snow layers observed during the SHEBA cam- 
paign (Sturm et al., 2002). Sturm et al. (1995) gave a rough 
generalization of snow kinds in different terrain types and 
suggested that the snow found at tundra is the closest substi- 
tute for what would be estimated to be found over floating 
ice. Accordingly, the values of snow density applied in this 
study are considered to be representative of the typical verti- 
cal profiles of snow density in and around the Arctic Ocean 
in spring. That is generally composed of wind slab, recent 
and new snow, fine-grained snow, wind slab, and depth hoar 
from top to bottom by turns (Sturm et al., 2002). In the es- 
timation of total-layer sBC, we presume that the uncertain- 


ties are entirely represented by the standard deviations of 
observed snow densities. We take the depth-weighted aver- 
age of the concentration values in each layer as the estimated 
vertically-integrated concentration when the value of surface 
sBC was close to or less than that of subsurface sBC. All 
surface and subsurface observations used in this study and 
the derived vertically-integrated values are listed in Table 1. 
In order to avoid influence from the surrounding urban envi- 
ronment, Table 1 only includes the estimated values of sBC 
in the Arctic Ocean and its adjacent coastal regions, exclud- 
ing ones far away from the Arctic Ocean. The vertically- 
integrated values larger than 70ngg _l taken in some sites, 
were also excluded because large difference between surface 
and subsurface concentrations was due to local contamina- 
tion, which may bring great uncertainties. The concentrations 
of sBC shown in Table 1 have been measured either by the 
“spectrophotometry” method (for most values in Doherty et 
al., 2010) or the “thermo-optical NIOSH 5040” method (for 
the values in this study and Forsstrom et al., 2009). Both of 
the methods discriminate BC and non-BC fractions of ab- 
sorption in the measuring process. It is generally thought that 
the “spectrophotometry” method is better in the albedo ef- 
fect analysis. However, in spatial distribution analysis, both 
of them could imply real variability because discrepancies 
among measured amounts greatly exceed the estimated er- 
rors induced by different measurement methods (Hansen et 
al., 2004). 
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Table 1 . Present-day sBC observations in and around the Arctic Ocean in spring. The estimated sBC in this study were calculated from 
surface and subsurface concentrations with the snow densities obtained in SHEBA field campaign. Most of the surface and subsurface values 
listed in this table are from Doherty et al. (2010) and others are from several field campaigns known to us. The measurements within a single 
model grid (“15-20”, “32-33”, “44-45”, “46-49”, “53-63”, “64-70”, “71-73”) are averaged before the comparison with modeled sBC. 


NO. 

Measurement 

Region 

Lat (N) 

Lon (E) 

Measured sBC 
in surface 
snow (ngg -1 ) 

Measured sBC 
in subsurface 
snow (ngg -1 ) 

Measured sBC 
in whole layer 
snow (ngg -1 ) 

Measurement 

period 

Reference 

Estimated 
sBC (ngg- 1 ) 

Uncertainty 
of Estimation 
(ngg -1 ) 

1 

Beaufort Sea 

73.000 

215.000 

9.94 

8.63 

- 

Apr 2007 

Doherty et al. (2010) 

8.89 

0.08 

2 

Beaufort Sea 

75.340 

224.343 

- 

- 

8.00 

Apr 2007 

APLIS/SEDNA campaign 

8.00 

- 

3 

Arctic Ocean 

82.880 

205.550 

- 

- 

7.00 

Apr-May 2006 

NPEO (Field campaign) 

7.00 

- 

4 

Arctic Ocean 

84.182 

47.507 

- 

- 

6.70 (n/a) 

May 2008 

Doherty et al. (2010) 

6.70 

- 

5 

Arctic Ocean 

84.363 

62.997 

- 

- 

10.90 (n/a) 

May 2008 

Doherty et al. (2010) 

10.90 

- 

6 

Arctic Ocean 

84.700 

296.520 

- 

- 

6.00 

Apr-May 2008 

Doherty et al. (2010) 

5.00 

- 

7 

Arctic Ocean 

85.897 

25.797 

- 

- 

5.80 (n/a) 

May 2008 

Doherty et al. (2010) 

5.80 

- 

8 

Arctic Ocean 

86.142 

45.342 

- 

- 

7.10 (n/a) 

Apr 2008 

Doherty et al. (2010) 

7.10 

- 

9 

Arctic Ocean 

86.300 

334.830 

- 

- 

6.75 

Apr-May 2008 

Switchyard campaign 

6.75 

- 

10 

Arctic Ocean 

88.133 

269.283 

- 

- 

4.60 (n/a) 

Spring 2007 

Doherty et al. (2010) 

4.60 

- 

11 

Arctic Ocean 

88.345 

3.407 

- 

- 

7.20 (n/a) 

Apr 2008 

Doherty et al. (2010) 

7.20 

- 

12 

Arctic Ocean 

88.560 

314.375 

- 

- 

5.00 

Apr-May 2008 

Switchyard campaign 

5.00 

- 

13 

Arctic Ocean 

88.850 

175.433 

- 

- 

9.20 (n/a) 

Apr 2006 

Doherty et al. (2010) 

9.20 

- 

14 

N. Pole 

89.150 

257.517 

- 

- 

6.00 (n/a) 

Spring 2007 

Doherty et al. (2010) 

6.00 

- 

15 

N. Pole 

89.300 

358.120 

- 

- 

4.85 (n/a) 

Spring 2008 

Doherty et al. (2010) 

4.85 

- 

16 

N. Pole 

89.400 

359.000 

- 

- 

2.33 (n/a) 

Spring 2008 

Doherty et al. (2010) 

2.33 

- 

17 

N. Pole 

89.500 

358.930 

- 

- 

2.51 (n/a) 

Spring 2008 

Doherty et al. (2010) 

2.51 

- 

18 

N. Pole 

89.700 

336.400 

- 

- 

3.00 

Apr-May 2007 

NPEO (Field campaign) 

3.00 

- 

19 

N. Pole 

89.900 

329.400 

- 

- 

4.20 (n/a) 

Spring 2008 

Doherty et al. (2010) 

4.20 

- 

20 

N. Pole 

90.000 

0.980 

- 

- 

4.00 (n/a) 

Spring 2008 

Doherty et al. (2010) 

4.00 

- 

21 

Canadian Arctic 

66.171 

255.626 

13.90 

6.68 

- 

Spring 2009 

Doherty et al. (2010) 

8.11 

0.46 

22 

Canadian Arctic 

67.878 

283.530 

9.40 

4.20 

- 

Spring 2009 

Doherty et al. (2010) 

5.23 

0.33 

23 

Canadian Arctic 

68.305 

255.913 

6.97 

10.30 

- 

Spring 2009 

Doherty et al. (2010) 

9.64 

0.21 

24 

Canadian Arctic 

68.568 

230.477 

9.10 

6.57 

- 

Spring 2009 

Doherty et al. (2010) 

7.07 

0.16 

25 

Canadian Arctic 

68.824 

264.711 

12.40 

7.30 

- 

Spring 2009 

Doherty et al. (2010) 

8.31 

0.32 

26 

Canadian Arctic 

68.986 

224.938 

7.15 

11.55 

- 

Spring 2009 

Doherty et al. (2010) 

10.68 

0.28 

27 

Canadian Arctic 

69.280 

282.954 

8.60 

5.13 

- 

Spring 2009 

Doherty et al. (2010) 

5.82 

0.22 

28 

Alaskan Arctic 

69.300 

216.200 

5.00 

- 

- 

Spring 2009 

Doherty et al. (2010) 

5.00 

- 

29 

Canadian Arctic 

69.635 

227.819 

10.00 

9.00 

- 

Spring 2009 

Doherty et al. (2010) 

9.20 

0.06 

30 

Canadian Arctic 

69.663 

250.904 

7.10 

5.70 

- 

Spring 2009 

Doherty et al. (2010) 

5.98 

0.09 

31 

Canadian Arctic 

69.895 

247.253 

6.00 

13.87 

- 

Spring 2009 

Doherty et al. (2010) 

12.31 

0.50 

32 

Canadian Arctic 

70.067 

235.027 

15.90 

10.50 

- 

Spring 2009 

Doherty et al. (2010) 

11.57 

0.34 

33 

Canadian Arctic 

70.067 

235.027 

12.10 

5.30 

- 

Spring 2009 

Doherty et al. (2010) 

6.65 

0.43 

34 

Canadian Arctic 

71.151 

280.752 

4.70 

3.43 

- 

Spring 2009 

Doherty et al. (2010) 

3.68 

0.08 

35 

Alaskan Arctic 

71.325 

203.567 

9.00 

- 

- 

Spring 2009 

Doherty et al. (2010) 

9.00 

- 

36 

Canadian Arctic 

72.341 

277.645 

9.50 

3.24 

- 

Spring 2009 

Doherty et al. (2010) 

4.48 

0.40 

37 

Canadian Arctic 

72.566 

259.193 

15.20 

6.30 

- 

Spring 2009 

Doherty et al. (2010) 

8.06 

0.57 

38 

Canadian Arctic 

75.497 

263.855 

9.97 

9.40 

- 

Spring 2009 

Doherty et al. (2010) 

9.51 

0.04 

39 

Canadian Arctic 

76.555 

255.268 

9.32 

8.65 

- 

Spring 2009 

Doherty et al. (2010) 

8.78 

0.04 

40 

Canadian Arctic 

76.633 

263.788 

11.70 

3.80 

- 

Spring 2009 

Doherty et al. (2010) 

5.37 

0.50 

41 

Canadian Arctic 

76.867 

274.786 

9.20 

5.83 

- 

Spring 2009 

Doherty et al. (2010) 

6.50 

0.21 

42 

Canadian Arctic 

80.083 

273.300 

12.00 

- 

- 

Spring 2007 

Doherty et al. (2010) 

12.00 

- 

43 

Western Russia 

67.631 

53.646 

19.00 

8.00 

- 

Spring 2007 

Doherty et al. (2010) 

10.18 

0.70 

44 

Western Russia 

73.381 

81.429 

17.60 

21.65 

- 

Spring 2007 

Doherty et al. (2010) 

20.85 

0.26 

45 

Western Russia 

73.428 

81.481 

12.00 

27.00 

- 

Spring 2007 

Doherty et al. (2010) 

24.03 

0.95 

46 

Western Russia 

72.176 

102.839 

83.60 

34.20 

- 

Spring 2007 

Doherty et al. (2010) 

44.00 

3.14 

47 

Western Russia 

72.211 

102.933 

23.20 

27.90 

- 

Spring 2007 

Doherty et al. (2010) 

26.97 

0.30 

48 

Western Russia 

72.277 

103.102 

45.80 

40.40 

- 

Spring 2007 

Doherty et al. (2010) 

41.47 

0.34 

49 

Western Russia 

72.244 

103.019 

62.10 

40.90 

- 

Spring 2007 

Doherty et al. (2010) 

45.10 

1.35 

50 

Eastern Russia 

72.054 

128.523 

87.90 

17.60 

- 

Spring 2008 

Doherty et al. (2010) 

31.54 

4.47 

51 

Eastern Russia 

74.065 

128.872 

13.00 

26.00 

- 

Spring 2008 

Doherty et al. (2010) 

23.42 

0.83 

52 

Eastern Russia 

71.649 

127.894 

82.00 

25.80 

- 

Spring 2008 

Doherty et al. (2010) 

36.94 

3.57 

53 

Eastern Russia 

68.631 

160.369 

83.20 

15.50 

- 

Spring 2008 

Doherty et al. (2010) 

28.92 

4.30 

54 

Eastern Russia 

68.649 

160.487 

80.00 

14.00 

- 

Spring 2008 

Doherty et al. (2010) 

27.09 

4.20 

55 

Eastern Russia 

68.663 

160.592 

88.30 

18.50 

- 

Spring 2008 

Doherty et al. (2010) 

32.34 

4.44 

56 

Eastern Russia 

69.043 

161.122 

55.60 

23.50 

- 

Spring 2008 

Doherty et al. (2010) 

29.87 

2.04 

57 

Eastern Russia 

69.044 

161.123 

53.15 

19.90 

- 

Spring 2008 

Doherty et al. (2010) 

26.49 

2.11 

58 

Eastern Russia 

69.032 

161.201 

50.00 

26.00 

- 

Spring 2008 

Doherty et al. (2010) 

30.76 

1.53 

59 

Eastern Russia 

69.019 

161.278 

39.30 

18.37 

- 

Spring 2008 

Doherty et al. (2010) 

22.52 

1.33 

60 

Eastern Russia 

69.020 

161.279 

52.40 

31.90 

- 

Spring 2008 

Doherty et al. (2010) 

35.97 

1.30 

61 

Eastern Russia 

68.737 

161.521 

57.70 

15.53 

- 

Spring 2008 

Doherty et al. (2010) 

23.89 

2.68 

62 

Eastern Russia 

68.719 

161.572 

53.00 

18.00 

- 

Spring 2008 

Doherty et al. (2010) 

24.94 

2.23 

63 

Eastern Russia 

68.700 

161.623 

42.05 

8.30 

- 

Spring 2008 

Doherty et al. (2010) 

14.99 

2.15 

64 

Eastern Russia 

69.869 

169.302 

11.00 

10.00 

- 

Spring 2008 

Doherty et al. (2010) 

10.20 

0.06 

65 

Eastern Russia 

69.782 

169.720 

11.90 

8.07 

- 

Spring 2008 

Doherty et al. (2010) 

8.83 

0.24 

66 

Eastern Russia 

68.930 

170.713 

14.30 

9.30 

- 

Spring 2008 

Doherty et al. (2010) 

10.29 

0.32 

67 

Eastern Russia 

69.330 

170.856 

11.50 

15.30 

- 

Spring 2008 

Doherty et al. (2010) 

14.55 

0.24 

68 

Eastern Russia 

69.119 

170.858 

14.00 

11.00 

- 

Spring 2008 

Doherty et al. (2010) 

11.59 

0.19 

69 

Eastern Russia 

69.022 

170.918 

23.20 

10.80 

- 

Spring 2008 

Doherty et al. (2010) 

13.26 

0.79 

70 

Eastern Russia 

69.195 

170.946 

13.20 

12.00 

- 

Spring 2008 

Doherty et al. (2010) 

12.24 

0.08 

71 

Eastern Russia 

69.571 

171.015 

9.50 

14.88 

- 

Spring 2008 

Doherty et al. (2010) 

13.81 

0.34 

72 

Eastern Russia 

69.524 

171.310 

13.00 

13.00 

- 

Spring 2008 

Doherty et al. (2010) 

13.00 

0.00 

73 

Eastern Russia 

69.478 

171.605 

18.00 

14.43 

- 

Spring 2008 

Doherty et al. (2010) 

15.14 

0.23 

74 

Western Svalbard 

78.910 

11.720 

- 

- 

3.20 

Mar-Apr 2007 

Forsstrom et al. (2009) 

3.20 

- 

75 

Western Svalbard 

78.870 

12.460 

- 

- 

1.70 

Apr 2007 

Forsstrom et al. (2009) 

1.70 

- 

76 

Western Svalbard 

79.000 

14.000 

- 

- 

1.40 

Apr 2007 

Forsstrom et al. (2009) 

1.40 

- 

77 

Eastern Svalbard 

77.897 

18.302 

- 

- 

9.80 

Mar 2007 

Forsstrom et al. (2009) 

9.80 

- 

78 

Eastern Svalbard 

78.750 

17.580 

- 

- 

6.60 

Mar 2007 

Forsstrom et al. (2009) 

6.60 

- 

79 

Eastern Svalbard 

79.910 

25.000 

- 

- 

6.50 

Apr 2007 

Forsstrom et al. (2009) 

6.50 

- 

80 

Northern Greenland 

81.000 

301.000 

- 

- 

2.60 (n/a) 

Spring 2008 

Doherty et al. (2010) 

2.60 

- 

81 

West of Greenland 

77.450 

299.500 

- 

- 

3.30 (n/a) 

Spring 2008 

Doherty et al. (2010) 

3.30 

- 

82 

West of Greenland 

76.40 

292.30 

- 

- 

4.30 (n/a) 

Spring 2008 

Doherty et al. (2010) 

4.30 

- 
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Table 2. Ratios of wet deposition to wet plus dry deposition in the GISS model from September to May 2007-2009. 



Sep 

Oct 

Nov 

Dec 

Jan 

Feb 

Mar 

Apr 

May 

North of 66° N 

0.95 

0.94 

0.92 

0.91 

0.90 

0.90 

0.92 

0.93 

0.92 

Northern Hemisphere 

0.78 

0.78 

0.77 

0.76 

0.74 

0.75 

0.76 

0.77 

0.76 

Worldwide 

0.81 

0.81 

0.81 

0.81 

0.81 

0.82 

0.82 

0.82 

0.81 


3 Model description 

The model GISS-PUCCINI is the NASA Goddard Institute 
for Space Studies (GISS) atmospheric composition and cli- 
mate model. It consists of the model for Physical Under- 
standing of Composition-Climate Interactions and Impacts 
(PUCCINI) (Shindell et al., 2006). which is fully embedded 
in the GISS modelE climate model (Schmidt et al., 2006). 
The atmospheric model version used here, GISS-E2, is that 
used for the Coupled Model Intercomparison Project (CMIP) 
phase 5 simulations in support of the IPCC fifth assessment 
report (AR5), and contains updates to the physics relevant 
to aerosols, including the ability to represent multiple down- 
drafts and updrafts in convective systems, while the black 
carbon model is unchanged from Koch and Hansen (2005). 
The model was run at 2° latitude by 2.5° longitude horizon- 
tal resolution with 40 vertical layers. Simulations were per- 
formed in a chemical transport model-like mode for 1995- 
2009 using observed sea surface temperatures (Rayner et al., 
2003) and linear relaxation of winds towards either NCEP 
or MERRA reanalysis (Kalnay et al., 1996; Rienecker et 
al., 201 1) to facilitate comparison with observations. Results 
are analyzed for 2006-2009. Repeating year 2000 monthly- 
varying emissions were used from the data set assembled for 
the AR5 simulations (Lamarque et al., 2010), with the excep- 
tion of biomass burning emissions which were monthly- and 
annually-varying emissions from the Global Fire Emission 
Database (GFED) version 3 (van der Werf et al., 2010). The 
emissions are respectively 3.5 Tg yr -1 from biomass burning 
and 5.3Tgyr -1 from fossil fuel and biofuel combustion. In 
this model, wet deposition depends upon solubility and on 
transport within convective plumes, scavenging within and 
below updrafts, rainout within both convective and large- 
scale clouds, washout below precipitating regions, evapora- 
tion of falling precipitation, and both detrainment and evap- 
oration from convective plumes. Wet deposition is the domi- 
nant removal process in this model. Its fraction of total depo- 
sition worldwide and over the Arctic has been shown in Ta- 
ble 2, from which we can see the values worldwide and over 
the Arctic regions are consistent with those given in other 
models of 78% to >95% (Textor et al., 2006; Huneeus et 
al., 201 1). Dry deposition is calculated using a resistance-in- 
series model coupled to a global, seasonally varying vegeta- 
tion dataset. More details are available in (Koch et al., 1999; 
Shindell et al., 2001). 


4 Initial field in the Arctic Ocean 

Another issue that needs to be considered is the initial snow 
cover (depth and density) and sBC distribution over multi- 
year sea ice at the beginning of the accumulation and de- 
position season (in this study - September). The initial field 
represents the snow and BC content within it that survived 
melt season. Because thermodynamic and dynamic mod- 
els of sea ice cover are limited in their ability to simu- 
late the time-varying ice conditions during the summer, we 
assume the initial snow cover to be that of the Septem- 
ber climatological snow conditions in Warren et al. (1999), 
which was also applied in the estimation of snow depth 
in Kwok Kwok and Cunningham (2008). The observations 
of snow and ice density (recent snow: 0.102 + 0.019 gem -3 
surface [<3cm]: 0.291 + 0.056 gem -3 subsurface [>3m]: 
0.333 + 0.02 gem -3 ) in the third Chinese Arctic expedition 
route (79.8-85.4° N, 144.1-170.1° W) in summer 2008 have 
been applied in the estimation of snow water equivalent 
(SWE) and in the calculation of vertically- integrated con- 
centrations of sBC from surface and subsurface measure- 
ments. In this study, we maintain a separate record of the ini- 
tial conditions at each grid cell. Table 3 summarizes the orig- 
inal observations and estimated values of sBC in the Arctic 
Ocean in summer 2005, 2008 and 2010. These measurements 
range from 74° N to 89° N, all of them are used to represent 
the climatological distribution of sBC over the multiyear sea 
ice in summer. The entire Arctic Ocean is divided into three 
latitudinal bands- south of 80° N, 80-85° N and the Arctic 
Ocean center, the distribution of sBC in each band is repre- 
sented by the mean values of the observations located within 
it. In this case, the initial fields are not time-varying and are 
identical for the results shown in this work. 

5 Intercomparison between model results and 
observations 

Present observations of sBC show sketchy but identifiable 
variation in spatial distribution, with maximum values in the 
Russian Arctic and much lower values in other regions. From 
the simulations, calculated from the model dry and wet de- 
position in spring 2007-2009, we can see that model val- 
ues are also higher over the Eurasian than the Canadian 
and Alaskan Arctic, followed by the Arctic Ocean, with 
the lowest in Greenland (Fig. 2). There is an apparent gra- 
dient from around the coast to the central Arctic Ocean, 
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Table 3. Present-day sBC observations in the Arctic Ocean in summer. The estimated sBC of whole layer at each measurement site is 
calculated from surface and subsurface BC concentrations with snow densities observed in the 3rd Chinese Arctic expedition in summer 
2008. Samples were collected on the sea ice surface, of melting sea ice, aged snow and newly-fallen snow from late July through late 
September. Thus, there are large uncertainties in current analysis, and present observations can only give a general and rough indication of 
the spatial distribution of sBC in summer. As the initial field of BC deposition, the Arctic Ocean is divided into several latitude zones in this 
study and each zone is characterized by one value averaged from the observations within this area. 


NO. 

Lat (N) 

Lon (E) 

Measured sBC 
in surface 
snow (ngg -1 ) 

Measured sBC 
in subsurface 
snow (ng g - 1 ) 

Measured sBC 
in whole layer 
snow (ng g - 1 ) 

Measurement 

period 

Reference 

Estimated 
sBC (ngg -1 ) 

Uncertainty of 

Estimation 

(ngg -1 ) 

1 

73.7 

192.59 

- 

- 

5.7 

Jul 2010 

This study 

5.7 

- 

2 

73.72 

192.72 

- 

- 

6.5 

Jul 2010 

This study 

6.5 

- 

3 

73.74 

193.04 

- 

- 

15.9 

Jul 2010 

This study 

15.9 

- 

4 

75 

200.01 

- 

- 

11.4 

Jul 2010 

This study 

11.4 

- 

5 

75 

202.5 

- 

- 

2.8 

Jul 2010 

This study 

2.8 

- 

6 

75.02 

199.98 

- 

- 

2.9 

Jul 2010 

This study 

2.9 

- 

7 

75.03 

200.52 

- 

- 

7.4 

Jul 2010 

This study 

7.4 

- 

8 

75.71 

222.8 



16 

Aug 2005 

Perovich et al. (2009) 

16 


9 

75.908 

219.411 

10.2 

10.2 

- 

Aug 2008 

Doherty et al. (2010) 

10.2 

0.0 

10 

75.908 

219.411 

12.2 

8.5 


Aug 2008 

Doherty et al. (2010) 

9.7 

0.2 

11 

75.99 

203.05 

- 

- 

41.6 

Jul 2010 

This study 

41.6 

- 

12 

76 

203.97 

- 

- 

10.2 

Jul 2010 

This study 

10.2 

- 

13 

76.81 

199.48 

- 

- 

40.2 

Aug 2010 

This study 

40.2 

- 

14 

77.04 

200.18 

- 

- 

22.1 

Aug 2010 

This study 

22.1 

- 

15 

77.04 

200.18 

- 

- 

9.5 

Aug 2010 

This study 

9.5 

- 

16 

77.04 

200.18 


- 

8.5 

Aug 2010 

This study 

8.5 

- 

17 

77.702 

213.398 

10 

- 

- 

Jul 2008 

Doherty et al. (2010) 

10.0 

- 

18 

77.9 

200.13 

- 

- 

39.2 

Aug 2010 

This study 

39.2 

- 

19 

77.98 

200.36 

- 

- 

29.8 

Aug 2010 

This study 

29.8 

- 

20 

77.98 

200.36 

- 

- 

10.8 

Aug 2010 

This study 

10.8 

- 

21 

78 

220.42 

- 


15 (n/a) 

Aug 2008 

Doherty et al. (2010) 

15.0 


22 

78.01 

209.818 

4.2 

15.2 

- 

Aug 2008 

Doherty et al. (2010) 

14.3 

0.6 

23 

78.291 

183.321 

12.3 

- 

- 

Aug 2005 

Doherty et al. (2010) 

12.3 

- 

24 

78.387 

206.549 

21 

22.1 

- 

Aug 2008 

Doherty et al. (2010) 

21.6 

0.1 

25 

78.392 

206.377 

16.4 

9.3 


Aug 2008 

Doherty et al. (2010) 

10.9 

0.4 

26 

78.6 

216.5 

- 

- 

26 

Aug 2005 

Perovich et al. (2009) 

26.0 

- 

27 

79.83 

331.21 

- 

- 

1.2 (n/a) 

Aug 2008 

Doherty et al. (2010) 

1.2 

- 

28 

79.88 

333.99 

- 

- 

4.8 (n/a) 

Aug 2008 

Doherty et al. (2010) 

4.8 

- 

29 

79.988 

209.713 

28.9 

- 

- 

Aug 2008 

Doherty et al. (2010) 

28.9 

- 

30 

80.003 

209.656 

49.4 


- 

Aug 2008 

Doherty et al. (2010) 

49.4 


31 

80.081 

209.792 

23.5 

13.7 

- 

Aug 2008 

Doherty et al. (2010) 

18.3 

0.6 

32 

81.226 

182.805 

3.5 

- 

- 

Aug 2005 

Doherty et al. (2010) 

3.5 

- 

33 

81.226 

182.805 

4.1 

- 

- 

Aug 2005 

Doherty et al. (2010) 

4.1 


34 

81.723 

209.035 

5.7 

8.6 

- 

Aug 2008 

Doherty et al. (2010) 

7.2 

0.2 

35 

81.926 

210.071 

14.1 

- 

- 

Aug 2008 

Doherty et al. (2010) 

14.1 

- 

36 

81.997 

219.943 

4.3 

17.7 

- 

Aug 2008 

Doherty et al. (2010) 

14.4 

0.8 

37 

83.087 

185.329 

3 

- 

- 

Aug 2005 

Doherty et al. (2010) 

3.0 

- 

38 

83.087 

185.329 

4.3 

- 

- 

Aug 2005 

Doherty et al. (2010) 

4.3 

- 

39 

83.299 

188.112 

1.8 

- 

- 

Aug 2005 

Doherty et al. (2010) 

1.8 

- 

40 

83.299 

188.112 

5.2 

- 

- 

Aug 2005 

Doherty et al. (2010) 

5.2 

- 

41 

83.955 

216.808 

4.2 

- 

- 

Aug 2005 

Doherty et al. (2010) 

4.2 

- 

42 

83.955 

216.808 

11.8 

- 

- 

Aug 2005 

Doherty et al. (2010) 

11.8 

- 

43 

84.171 

209.005 

2.6 

- 

- 

Aug 2005 

Doherty et al. (2010) 

2.6 

- 

44 

84.171 

209.005 

7 

- 

- 

Aug 2005 

Doherty et al. (2010) 

7.0 

- 

45 

84.307 

210.918 

4.6 

- 

- 

Aug 2005 

Doherty et al. (2010) 

4.6 

- 

46 

84.309 

199.352 

2 

- 

- 

Aug 2005 

Doherty et al. (2010) 

2.0 

- 

47 

84.309 

199.352 

14.4 

- 

- 

Aug 2005 

Doherty et al. (2010) 

14.4 

- 

48 

84.311 

199.581 

2.2 

- 

- 

Aug 2005 

Doherty et al. (2010) 

2.2 

- 

49 

84.311 

199.581 

15.7 

- 

- 

Aug 2005 

Doherty et al. (2010) 

15.7 

- 

50 

85.122 

205.2 

18.8 


- 

Sep 2005 

Doherty et al. (2010) 

18.8 

- 

51 

85.938 

48.335 

6.4 

2.8 

- 

Sep 2005 

Doherty et al. (2010) 

4.5 

0.2 

52 

85.938 

48.335 

9.6 

- 

- 

Sep 2005 

Doherty et al. (2010) 

9.6 

- 

53 

86.657 

55.618 

3.7 

6.1 

- 

Sep 2005 

Doherty et al. (2010) 

5.0 

0.1 

54 

87.472 

57.588 

4.2 

9.6 

- 

Sep 2005 

Doherty et al. (2010) 

7.1 

0.3 

55 

87.62 

155.876 

2.7 

- 

- 

Sep 2005 

Doherty et al. (2010) 

2.7 

- 

56 

87.62 

155.876 

30 

- 

- 

Sep 2005 

Doherty et al. (2010) 

30.0 

- 

57 

87.66 

150.902 

8 

- 

- 

Aug 2005 

Perovich et al. (2009) 

8.0 

- 

58 

88.056 

58.748 

2.3 

2.1 

- 

Sep 2005 

Doherty et al. (2010) 

2.2 

0.0 

59 

88.456 

146.532 

5.3 

- 

- 

Sep 2005 

Doherty et al. (2010) 

5.3 

- 

60 

88.456 

146.532 

22.4 

- 

- 

Sep 2005 

Doherty et al. (2010) 

22.4 

- 

61 

88.46 

213.47 

7 

- 

- 

Aug 2005 

Perovich et al. (2009) 

7.0 

- 

62 

88.813 

164.136 

18.2 

- 

- 

Sep 2005 

Doherty et al. (2010) 

18.2 

- 

63 

88.813 

164.136 

1.7 


- 

Sep 2005 

Doherty et al. (2010) 

1.7 

- 

64 

89.374 

270.912 

2.7 

13.3 

- 

Sep 2005 

Doherty et al. (2010) 

8.4 

0.6 

65 

89.482 

169.798 

7.1 

2.3 

- 

Sep 2005 

Doherty et al. (2010) 

4.5 

0.3 
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2007 


2008 


2009 



Fig. 2. Simulated sBC in and around the Arctic Ocean in spring 2007-2009. Black circles indicate the measurement sites for the correspond- 
ing years. 


especially in the Russian Arctic, but relatively little vari- 
ability over the Arctic Ocean. The mean concentration of 
sBC over the Arctic Ocean shows large interannual variabil- 
ity in spring from 2007 to 2009 (7.8-13.4 ngg -1 ). Point to 
point comparisons between measured sBC and model sim- 
ulations from NCEP and MERRA runs have been carried 
out respectively in the Arctic Ocean and surrounding seas, 
the Canadian and Alaskan Arctic, the Russian Arctic, Sval- 
bard and Greenland (Fig. 3). Model results have been inter- 
polated to each measurement site in the corresponding month 
and year of observations and are compared with the mean 
observed value in cases where more than one measurement 
is available within a 2 x 2.5 degree grid box. From Fig. 4, 
it can be seen that both of the simulated values from two 
model runs are comparable with most of present observa- 
tions over each Arctic sub-region. Rough agreement between 
the model and observations is found in the Arctic Ocean and 
Canadian and Alaskan Arctic sector, respectively with mean 
values of 8.2 + 0.9ngg _1 and 7.2+ 1.3ngg _1 in the model 
and 7.4 + 2.3 ng g ~ 1 and 7.8 + 2.4 ng g“ 1 in the observations. 
The concentrations of modeled sBC in Western Greenland 
in spring 2008 are also close to the observations; mean val- 
ues are 3.8 + 1.8 ng g -1 in the model and 3.8 + 0.7 ngg -1 in 
the observations. However, there is an apparent difference in 
Northern Greenland, which conforms to the abnormal peak- 
value in the central inland areas of Greenland (Fig. 3). More 
measurements would be needed to determine the actual dis- 
tribution of sBC there. Additionally, the values of modeled 
sBC are significantly smaller than observations in the Rus- 
sian Arctic sector in spring 2007. Previous studies indicate 
that anthropogenic influence was dominant in western Rus- 
sia in spring 2007, and open fire influence was dominant in 
Eastern Russia in spring 2008 (Doherty et al., 2010; Wang et 
al., 2011). Indeed, the modeled values in Eastern Russia are 
substantially higher in 2008 than in 2007, perhaps owing to 
the interannually varying biomass burning emissions. Hence 


an underestimation of biomass burning emissions in Russia 
may be the main reason of this simulation bias. 

The simulation values in Svalbard are within the range 
of observations and near to the median value. Analysis of 
previous measurements (81 samples) in Svalbard indicates 
that there exist obvious spatial variations, the concentrations 
of sBC in eastern Svalbard are significantly higher than the 
western side (Forsstrom et al., 2009) and the emissions from 
the western European Arctic may be dominant in this dis- 
tribution pattern. This spatial variation has also been repro- 
duced by the model even though the modeled concentrations 
in western Svalbard are a little higher than the observed con- 
centrations. 

Further study on the simulation from two model runs 
shows small difference over each Arctic sector. NCEP run 
performs slightly better in the Arctic Ocean, Svalbard and 
Greenland, MERRA run performs slightly better in Russia, 
Canadian and Alaskan Arctic, although they both show a sig- 
nificant underestimation in the Russian Arctic, and overesti- 
mation in Svalbard and Greenland (Table 4). The difference 
between the two model runs was far less than current simula- 
tion bias. Therefore, we suggest that the discrepancy between 
model results and observations comes from other sources and 
not from different meteorology fields applied. 

We also estimated the decrease in snow and sea ice albedo 
and resulting radiative forcing based on the sBC distribu- 
tion simulated by the current GISS-PUCCINI model. Fig- 
ure 4 shows the albedo reduction in and around the Arctic 
Ocean in spring 2007-2009 due to BC deposition. We es- 
timated the impact of sBC on snow and ice albedo accord- 
ing to the “spectrally averaged” albedo change by BC con- 
tent illustrated in Fig. 2 in Warren and Wiscombe (1985) 
and assume the snow grain radius is a constant of 100 pm 
(McConnell et al., 2007) with no significant aging. The re- 
sulting albedo reduction presents significant space-time vari- 
ations, with highest mean value of 1 .25 % in the Russian 
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Fig. 3. Histogram of modeled versus observed sBC for the Arctic Ocean and surrounding regions in spring 2007-2009. Model results are 
sampled for the month and year of observations and interpolated to corresponding measurement sites by bilinear interpolation methodology. 
Simulated values from NCEP and MERRA runs are shown both. 


2007 2008 2009 



Fig. 4. Modeled change in snow and ice albedo due to sBC in and around the Arctic Ocean (north to 66° N) in spring 2007-2009. 


Arctic, which was much larger than those in other Arctic re- 
gions ranging 0.39 % to 0.64%. The averaged value over the 
Arctic north of 66° N is 0.4-0. 6 % from spring 2007-2009, 
lower than the estimates of Park et al. (2005), but compara- 
ble with Jacobson et al. (2004). On the basis of this result 
of albedo decrease, we evaluate the radiative forcing caused 
by sBC with the NCEP downwelling surface solar radiation. 
Result shows that the radiative forcings from BC deposition 
to the Arctic snow and ice (north to 66° N) are 0.7 W m -2 , 

1.1 W m~ 2 and l.OWm -2 , respectively in spring 2007, 2008 
and 2009. Wang et al. (2011) reported a radiative forcing of 

1 .2 W rrT 2 for spring 2007-2009 in the Arctic north to 60° N, 
similar to Flanner et al. (2007) and slightly larger than our 
values. 


Finally, we review the research on the potential source 
type/region of Arctic BC in previous studies. Rahn and Mc- 
Caffrey (1980) reported that the Arctic atmosphere is hazy 
in winter and spring, and indicated this phenomenon may 
be caused by fossil fuel burning, industrial and agricultural 
processes, by long-range transport of mid-latitude pollution 
products. Generally, Europe, North America and South Asia 
are considered to be the main contributors to present Arc- 
tic sBC, while contributions from biomass burning sources 
are also significant (Hegg et al., 2009). Matsui et al. (2011) 
suggests that biomass burning and anthropogenic sources in 
high-latitude were most important for the Arctic BC both 
in spring and summer, because most of BC from lower- 
latitude regions has been removed by wet deposition before 
arriving at the Arctic. Among these source regions, Europe 
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Table 4. Observation-to-model ratios between observations and model results from NCEP and MERRA runs over each Arctic sector. 



Arctic Ocean 

Canada and Alaska 

Russia 

Svalbard 

Greenland 

NCEP 

0.92 

1.09 

1.81 

0.84 

0.81 

MERRA 

0.89 

1.05 

1.77 

0.73 

0.78 


contributes more than 50 % of the Arctic BC loading mainly 
through transport within the lower troposphere in all seasons 
(Huang et al., 2010b), which is the most likely cause of the 
highest concentrations in the Russian Arctic and surrounding 
Seas. Eleftheriadis et al. (2009) found some differences in 
the seasonal contribution to the BC measured at Svalbard by 
Eurasian sources with the BC peak values mainly during the 
winter arriving from Siberian and Western Russian sources, 
while during the warm season North Europe including West 
Russia was the main contributing area. Hegg et al. (2009) 
also indicates that local emissions from fossil-fuel combus- 
tion make a significant contribution to the springtime high- 
latitude Arctic Ocean and some locations in Western Russia. 
North America contributes about 10-20% to the Arctic tro- 
posphere with least variations in the contribution from differ- 
ent vertical levels (Huang et al., 2010b; Shindell et al., 2008), 
and is also regarded as a main biomass source to Green- 
land, the Alaskan and the Canadian Arctic (McConnell et 
al., 2007; Hegg et al., 2009). Because of its remoteness from 
the main source regions, the concentrations in Greenland and 
the central Arctic Ocean are much lower than other Arc- 
tic regions. South Asia is also considered to be a signifi- 
cant contributor in the Arctic upper troposphere/lower strato- 
sphere during the springtime but it does not appear to con- 
tribute significantly to the deposited BC (Koch and Hansen, 
2005; Hirdman et al., 2010). In addition to the pollution from 
long range transport sources, the emissions from ships and 
cross-polar aircraft flight are also very important. Corbett 
et al. (2007) indicates that the impacts from local shipping 
are comparable with long range transport from lower latitude 
emission sources. Emission form international shipping can 
be up to 71000-160000 metric tons annually, representing 
about 2 % of global BC from all sources. It is also shown 
that 1.2GgBC has been emitted from shipping within the 
Arctic in 2004, most of which was derived from transit ves- 
sels, container and general cargo ships, and the rest from fish- 
ing vessels (Corbett et al., 2010a). The vessel activities are 
mainly concentrated in summer and autumn months, and are 
increasing dramatically with the rapid decrease of the Arc- 
tic sea ice. In this background, many approaches have been 
recommended to reduce PM emissions from diesel engines 
and ships, such as switching to low-sulfur fuels and engine 
process modifications (Corbett et al., 2010b). Cross-polar air- 
craft flight is also an important contributor to the Arctic BC. 
Whitt et al. (2011) point out that 2.35 Tg fuel were burned in 
the Arctic Circle (66.56° N —90° N) in 2006 and more than 
half of commercial aviation emissions occur in relatively sta- 


ble regions of the atmosphere and nearly one quarter occur 
in the stratosphere. In view of the validation against present 
spatially extensive observations of sBC, we suggest that the 
latest NASA GISS-E2-PUCCINI model would be a rational 
choice to further identify source regions of Arctic BC. 

6 Discussion and conclusions 

There are large uncertainties in current model simulations 
of BC in the Arctic. Various models give very different re- 
gional source contributions and different evaluation of radia- 
tive forcing due to BC deposition. It is useful to characterize 
model biases by comparing results with present field obser- 
vations, in order to reduce these uncertainties. In this study, 
spatially extensive comparisons between model results and 
present sBC observations are presented, to validate the lat- 
est GISS model deposition of BC in the Arctic regions. It 
is demonstrated that the latest GISS model performs well in 
large scale simulations of BC deposition, especially in the 
Arctic Ocean and Canadian Arctic sector. However, there is 
still apparent underestimation for both runs in a few measure- 
ment sites in Russian Arctic. Further comparison of modeled 
aBC with surface long-term observations at Barrow, Alert, 
Zeppelin and Nord stations shows that the latest GISS model 
has been greatly improved in simulating the seasonal vari- 
ations of the Arctic BC (Fig. 5), in comparison to previous 
version mentioned in Shindell et al. (2008), but still has sig- 
nificant underestimation in winter and spring, especially at 
Barrow station, where the bias can be up to one order of 
magnitude. This is in agreement with the results from ear- 
lier studies (Koch et al., 2009b; Schwarz et al., 2010), in 
which comparisons of vertical profiles of simulated BC with 
corresponding observations had been done over the Arctic 
Ocean, North American, Canada and Greenland in spring and 
during a flight line from Anchorage to the north pole (60— 
80° N) in winter. The results indicated that the simulated ver- 
tical profiles were comparable with observations in Green- 
land, Canada and Arctic Ocean in spring. However, there 
was significant bias in Alaska and nearby regions whether in 
spring or winter, with the model underestimating BC concen- 
trations in the middle and lower troposphere. It is generally 
thought that several factors may be the cause of this bias. 
Earlier studies pointed out that the lack of seasonal varia- 
tions from current simulations may be due to poorly sim- 
ulated deposition processes rather than transport processes 
(Schwarz et al., 2010; Huang et al., 2011; Quinn et al., 201 1). 
Liu et al. (2011) found that winter concentrations of BC 
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Fig. 5. Comparison of modeled monthly aBC concentrations in near surface layer versus measured concentrations for Alert (83° N, 62° W), 
Barrow (71°N, 157° W), Zeppelin (79° N, 12° E) and Nord (81° N, 16° W) station from January 2006 to December 2009. 


in the Arctic could increase throughout much of the tropo- 
spheric column by improving parameterization of BC aging 
and deposition processesin the model. Huang et al. (2011) 
suggested that the dry deposition process has little effect on 
the seasonal pattern of BC in the Arctic lower troposphere, 
the observed seasonality of BC in the Arctic troposphere is 
mainly due to the seasonal changes in aerosol removal by 
wet scavenging and seasonal injections of BC from Europe 
and the former USSR. Koch et al. (2009) also suggested that 
current simulation bias in Arctic BC may be due to lack of 
sufficient removal by precipitating clouds and low-level pole- 
ward transport. Therefore, we expect current model output 
could perhaps be enhanced by improving the physical param- 
eterization of wet deposition processes and biomass burning 
emissions in the Arctic. Of course, this still needs to be veri- 
fied by a series of sensitivity tests with the model. However, 
present BC observations, obtained by means of various mea- 
surement technologies, also have large uncertainties (Koch 
et al., 2009). In addition, aircraft measurements were actu- 
ally carried out during short flight campaigns, which basi- 
cally provide a “snapshot” of vertical profiles of BC concen- 
trations. This makes the comparison with simulated monthly 
BC distribution quite challenging since one emission event 
from a forest or grass fire in the region, could modify the 
vertical BC profile and greatly affect the results of aircraft 
measurements. 

It is especially worth noting that the largest simulation bi- 
ases of sBC appeared in the Russian Arctic (see Fig. 3). Al- 
though there are a few measurement sites included in the 
comparison from this region, we expect that biomass burn- 


ing emissions may be the main cause of the poor perfor- 
mance here. On the other hand, Arctic industrial emissions 
could also play a role and could have changed over the 2000s 
decade. There are two pieces of evidence to support this as- 
sumption: (1) biomass burning and anthropogenic sources in 
high-latitude were regarded as the most important contribu- 
tor to the Arctic BC (Matsui et al., 2011), and emission from 
biomass burning is dominant in the contribution to Arctic 
sBC (Hegg et al„ 2009). (2) Strong cooling of the air over 
the snow cover in Eurasia allows polluted air from northern 
Eurasia to penetrate the entire Arctic at low altitudes (Stohl 
et al., 2006), and Northern Eurasia is most likely to be the 
source region of the aBC in the lower troposphere in the 
Arctic in winter and early spring (Quinn et al., 2011). Fi- 
nal conclusions cannot be reached from the above analyses. 
Further model simulations of BC are needed to determine if 
emissions are responsible, and if so to show which type of 
emissions from what regions and which type of sources lead 
to the current bias found in the Russian Arctic. For example, 
running the model in a period when Russian forest fires were 
minimal could provide some insight. In addition to the inher- 
ent problems in the model, the approximations and hypothe- 
ses used in the calculation of vertically-integrated concentra- 
tion of sBC may be another reason for the underestimation. 
We applied climatological observations of snow density ob- 
tained in the SHEBA campaign to simulate the actual values 
of snow density in different depths at each measurement site 
and this may lead to much larger uncertainties in the results 
obtained for the Russian Arctic compared with other Arc- 
tic regions, as significant variations in the vertical profiles 
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are observed there. More measurements of snow density in 
various snow depths are needed to reduce this uncertainty. 
Finally, more measurements with large spatial coverage in 
the Arctic would also be particularly useful for model im- 
provement. For instance, more aircraft observations for BC 
are needed over Eurasia, the oceans and the biomass burning 
regions. 
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